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(Dated: February 1, 2008) 

We calculate cross sections for low energy elastic exciton-exciton scattering within the effec- 
tive mass approximation. Unlike previous theoretical approaches, we give a complete, non- 
perturbative treatment of the four-particle scattering problem. Diffusion Monte Carlo is used 
to calculate the essentially exact energies of scattering states, from which phase shifts are 
determined. For the case of equal-mass electrons and holes, which is equivalent to positronium- 
positronium scattering, we find a s = 2.10^ for scattering of singlet-excitons and a s = 1.5a x 
for triplet-excitons, where a x is the excitonic radius. The spin dependence of the cross sec- 
tions arises from the spatial exchange symmetry of the scattering wavefunctions. A significant 
triplet-triplet to singlet-singlet scattering process is found, which is similar to reported effects in 
recent experiments and theory for excitons in quantum wells. We also show that the scattering 
length can change sign and diverge for some values of the mass ratio mh/m e , an effect not seen 
in previous perturbative treatments. 

PACS numbers: 07.05. Tp; 71.35.Cc; 02.70.Lq 



I. INTRODUCTION 

Excitons in semiconductors have been the subject 
of many experimental and theoretical investigations of 
Bose condensation. Low-energy exciton-exciton interac- 
tions are characterized by the exciton-exciton scattering 
length, a s , which determines the thermodynamics of a 
low density gas and is crucial for modeling the thermal- 
ization time of a dilute exciton gas. Despite its impor- 
tance, the exciton-exciton scattering length is an elusive 
quantity, being difficult to measure experimentally or to 
estimate theoretically. 

As is well-known in atomic physics, scattering lengths 
can be extremely sensitive to the details of the interac- 
tions between particles. In particular, the existence of 
a weakly bound or nearly bound state causes the scat- 
ter length to become quite large. Therefore, a priori one 
should suspect that exciton-exciton scattering may be a 
very material dependent property of semiconductors. Re- 
liable theoretical predictions of exciton-exciton scatting 
lengths require both a very accurate Hamiltonian for the 
semiconductor, and an accurate solution to the (four- 
particle) scattering problem. In this paper we provide 
an essentially exact solution to exciton-exciton scattering 
for a commonly used single-band effective mass Hamilto- 
nian. This solution allows us to study three important 
questions: (1) how sensitive is the scattering length to the 
mass ratio m e /mh, (2) how does the scattering length de- 
pend on spin states (singlet or triplet) of the scattering 
excitons, and (3) to what degree can inter-exciton ex- 
change of electrons or holes cause excitons to scatter into 
different spin states? This calculations also serve as a 



benchmark for the single band limit of more complicated 
scattering Hamiltonians. 

One experimental method for measuring the exciton 
scattering cross section is to look at line width broad- 
ening of the recombination spectra in a gas of excitons. 
Collisions between excitons increase the line width, caus- 
ing the line width to depend on the exciton-exciton scat- 
tering rate, nav, where n is the density and v is a typical 
exciton velocity. Extracting cross sections from a line 
width requires that (I) the density and velocity distribu- 
tion are known, and (2) elastic scattering is the fastest 
process. As discussed below, CU2O is a good material 
for corctparison to the model studied in this work. Snoke 
et a/lira have performed such experiments on CU2O and 
have found a line width broadening that suggests an up- 
per bound of Aa x on the scattering length. Although our 
simulations do not exactly model CU2O, we will compare 
our results to this value. 



II. THEORETICAL BACKGROUND 

Theoretical approaches to this problem start with the 
effective mass approximation, in which the system under 
consideration consists of two electrons, labeled 1 and 2, 
and two holes, labeled a and b. The Hamiltonian is 



H = -AiV? - - \ a V 2 a - X b V( 
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where A = Ti 2 /2m. The Hamiltonian has symmetry un- 
der exchange of electrons and exchange of holes, so eigen- 
states may be denoted by two exchange quantum num- 
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bers. The s-wave states are symmetric under exchange of 
excitons; a condition which is satisfied by states <p ++ and 
</>~~, where the +(— ) signs refer to (anti)symmetry under 
exchange of electrons and holes, respectively. Although 
this Hamiltonian is a well-accepted model for exciton- 
exciton scattering, we should point out a few of its defi- 
ciencies. For small excitons, such as those in CU2O, that 
have radii not much larger than the lattice spacing, non- 
parabolic terms in the kinetic energy and other correc- 
tions to the potential energy may be necessary. For many 
semiconductors, such as Si and Ge, the valence band is 
a mixture of three bands and cannot be described by 
a single parabolic band. In the case of Cu 2 0, the va- 
lence band is the parabolic spin-orbit split off band, and 
there are fewer complications. Interband exchange (vir- 
tual electron-hole recombination) is an important effect 
that has been neglected, and could be modeled by an 
additional spin-dependent potential term. 

This Hamiltonian also describes a family of scatter- 
ing processes for other particles, including hydrogen- 
hydrogen, positronium-positronium, and muonium- 
muonium scattering. The equal mass case is at 
an extreme (positronium scattering), where the Born- 
Oppenheimer approximation is the least applicable. 

There have been several theoretical estimates of 
exciton-exciton-,|-Kcattering for bulk systems&Qtil and 
quantum wellstffi as n well as calculations on biexciton- 
biexciton scattering,!! Only the bulk, elastic scattering 
calculationsSu are directly comparable to the results of 
this paper, but the techniques presented here could be 
generalized to the other scattering problems. Also, the 
results presented here provide a benchmark for evaluat- 
ing the approximations used in other theoretical treat- 
ments, and could lend insight into the reliability of the 
approximations in more complicated situations. 

One standard theoretical approach is diagrammatic 
perturbation theory, as presented in the-work of Keldysh 
and Kolsova and Haug and Hanamura.E3 They estimate 
the exciton-exciton scattering matrix as arising from a 
single term, (ki +q 7 k 2 -q\H- mt \ki, k 2 ), where \ki,k 2 ) rep- 
resents a state of two noninteracting excitons with mo- 
mentum k\ and k 2 and H luij is the inter-exciton Coulomb 
interaction. This method gives an estimate of a s = ya x 
(independent of the mass ratio), where a x = m~ x + m^ 1 
is the exciton radius, but it is an uncontrolled approxi- 
mation which may have limited validity in the low energy 
limit. One serious drawback of the method is that it does 
not include effects of the biexciton in the scattering. As 
we show later, biexciton vibrational states cause strong 
dependence of the scattering length on the mass ratio 
m e /mh, which is not captured by the low order pertur- 
bation theory. 

A second common approach was developed by Elko- 
moss and Munclryta and uses an effective exciton-exciton 
potential defined by V c ff(R) = (0f (R)\H \cj)f (R)), where 
4>f(R) is the wavefunction for two free excitons a dis- 
tance R apart. The effective potential V e s arises from the 



Hartree term and is used in a two-particle central-field 
calculation. While an exciton-exciton scattering pseudo- 
potential would be a very useful tool, this approximate 
form has some serious drawbacks. Among its deficiencies 
are a lack of correlation, no van der Waals attraction, 
a failure to reproduce biexciton states, and a vanishing 
interaction potential for m e = m^. The cross sections 
calculated by this method are small and lack qualitative 
agreement with the results of the present work. 

Some insight into exciton-exciton scattering can be 
gained by considering the bound states, biexcitons. Since 
the number of bound states, enters in the phase 
shift at zero energy, 5(0) = ttN-q, it is necessary that 
a good computation method for low energy scattering 
be able accurately calculate biexciton binding energies.til 
For the mass ratios considered (and far beyond, includ- 
ing deuterium) the biexcitons cannot bind in the <p~" 
states, so biexcitons in rotational s states always have 
<f> ++ wavefunctions. Detailed theoretical descriptions of 
biexcitons can be found in Ref. [jl2| . The equal mass 
case was shown to have a bound biexciton by Hyllerass 
and Ore using a variational argument ,113 and a better 
variational estimate of theJpinding energy was given by 
Brinkman, Rice and Bell,ta who found Eb = 0.029^, 
where E x = 0.5/(m f r 1 + m^ 1 ) is the exciton binding 
energy. However, because of the importance of corre- 
lation energy, the latter variational treatment was miss- 
ing half of the biexciton binding energy, as s h own by 
diffusion Monte Carlo (DMC) calculations,!!^ which 
find E B = 0.06404(4)^. DMC is a quantum Monte 
Carlo (QMC) method that uses a random walk to project 
out the ground state wavefunction from a variational 
wavefunction, in order to stochastically sample the exact 
ground state energy. The success of DMC for calculating 
biexciton energies has been a motivation for its use in the 
present scattering calculations. 

III. METHOD: QUANTUM MONTE CARLO 
CALCULATION OF SCATTERING 

A R-Matrix approach and scattering boundary 
conditions 

The i?-Matrix approach to scattering is to examine 
the standing waves of the system.._As shown by Carl- 
son, Papdharipande, and WiringaEZl and Alhassid and 
Koonin,Ea by fixing nodes in the standing waves the scat- 
tering problem may be cast as a stationary state problem 
suitable for QMC methods. For an elastic scattering pro- 
cess, we label the distance between the products by R, 
and the reduced mass of the products by m r . In exciton- 
exciton scattering there is a subtlety in the definition of 
R due to inter-exciton exchange, which we will address 
below in our discussion of the exciton-exciton scattering 
wavefunctions. Nonetheless, for large separation R, the 
relative motion of the products is free-particle like, so the 
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many-body wavefunction depends on R as 



1. 



4> oc sin[fci? — -In + Si(k)], 



(2) 



where I is the relative angular momentum, k is the scat- 
tering momentum, and <5;(fc) is the phase shift. If we 
constrain the wavefunction to have a node at a large ex- 
citon separation i?„, we find a discrete energy spectrum 
E a (R n ) 1 which may be computed by ground state or ex- 
cited state methods, such as DMC. Each choice of R n 
gives a spectrum of states <fi a , a — 1, 2, . . . , with energies 
E a that determine values of Si(k), 



Si{k a ) 



foot Rn 



1 



lir 



(3) 



where k a = \/2m r E a . The scattering matrix elements 
are determined by the phase shifts, 



Si(k) 



1 

2ik 



exp[2i<S;(fc)] - 1. 



(4) 



Carlson has also proposed fixing the logarithmic 
derivative of the wavefunction At the boundary instead 
of setting wavefunction to zero, 113 



S- V* 



(5) 



R=R b 



were S is the normal to the boundary surface, at a 
fixed radius Rb, and /3 parameterizes the boundary con- 
dition. This formulation has the advantage of separat- 
ing the choice of simulation size Rb (subject to Rb lying 
in the asymptotic region), from the sampling of energy, 
which is handled by varying /3, and is particularly well 
suited for finding the scattering length. The application 
to VMC calculations is straight-forward, but preserving 
the boundary condition in DMC calculations requires a 
method of images. E2l The results presented here do not 
use the logarithmic-derivative boundary condition. 



where dk(t) is the fcth eigenvector with eigenvalue Ak(t) 
and the matrices N and H are the overlap and Hamilto- 
nian matrices in our trial basis, given by 

N af3 (t) = J dR 1 dR 2 r a (R 2 )e- tH f p {R 1 ), (7) 
H a p(t) = J dR 1 dR 2 H f* a {R 2 )e- tH f {Ri). (8) 

The parameter t is the projection time. The eigenvalues 
Afc (t) are energy eigenvalues E% within the Hilbert space 
spanned by the projected trial functions {e~ tH / 2 f a }, and 
approach the exact energy eigenvalues in the limit of large 
t. 

The matrices N and H are sampled with random 
walks, using a guiding function ip which must be positive 
everywhere. The guiding function must have significant 
overlap with all basis functions, and should be optimized 
to decrease the variance of the sampled matrices. At each 
step i + 1 in the random walk, the coordinates Ri+i of 
the particles are updated using 

R l+1 =R t + rA^~ 1 V^(i? l ) + (2rA)2 Xl , (9) 

where Xi is a normally distributed random variate with 
zero mean and unit variance and r is the time step. In the 
limit of small t, Eq. @ describes a process for sampling 

The matrix element of e tH is estimated by integrat- 
ing the local energy of the guiding function, Eljp(R) = 
i/j^ 1 (R)Hip(R), along the random walk, 

-t + Y. 1 i[£^(R 3 -)+i^*(A,-+i)] 
W n , n+k =e i=» . (10) 

The estimators for matrices N and H are 
1 P 

n a p{kT) = - J2K(^)W t , l+k F l3 (R l+k ) (11) 
1 p 

h a p{kr) = -22K(Ri)Wi,i+kFp(Ri +k )E L0 {Ri + ^l2) 



B Calculation of excited states 

The use of excited states 4> a {a > 1) is necessary when 
there is a bound state and, more generally, when the scat- 
tering state being studied has its first node before the 
asymptotic region is reached. We use a QMC method 
to calawlate excited states developed by Ceperley and 
Bernu,Eil to adapt VMC and DMC methods for a Hilbert 
space of several low energy wavef unctions. A set of m 
trial wavefunctions is chosen, f a , a = 1, . . . 7[ m. The gen- 
eralized eigenvalue equation to be solved is,^ 



Ei 

0=1 



H a p(t) - A k (t)N a0 ]d k p(t) = 0, 



(6) 



where F a = f a /ip and E Lf3 = f f3 1 (R)Hf /3 (R) are the 
local energies of the trial basis states. 



C Form of exciton-exciton scattering wavefunctions 

We now discuss the form for the exciton scattering 
functions and j^ . As mentioned before, inter- 

exciton exchange of particles complicates the definition of 
exciton-exciton separation. There are two configurations 
for well-separated excitons, as shown in Fig. |l|. Configu- 
ration I has the electrons and hole paired as la, 26; and 
Configuration II as 16, 2a. We choose wavefunctions /* 
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Configuration I Configuration II 
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«<±j 



Q 2 1 Q r 1Ka . Q 



1a-2b _Y Y....!. 1 b.-2a 

0b b© £>a 



(|> ++ = ^ + (i)" (Symmetric s-wave) 
<|r~ = (I)' - (I)" (Anti-symmetric s-wave) 



FIG. 1: The two contributions to the exciton scatter- 
ing wavefunctions 4>, at large exciton-exciton separation. 
The asymptotic form of the S-wave scattering states are 
symmetric (4> ++ ) or antisymmetric (</> ) combinations 
of these configurations. This symmetric/ antisymmetric 
form is used for the trial wavefunctions / ++ and / , 
using / : and f u from Eq. (fTi]). 

and to represent these states, 



■I a 



(1: 



C 



- l + d q 



" T+d 



h r lb 



e-^e-^U a (r lb _ 2a ) 



f!2 



l + d„ 



c h r la 



a ^F, (13) 

where j,Cf,df,c g ,dg,Ch,dh, and parameters in the func- 
tion U a (r) are variational. These wavefunctions repre- 
sent two excitons in a relative s-wave state. Since these 
are not eigenstates of the exchange operator for elec- 
trons P e or holes, J\, we take linear combinations of 
the two for our trial wavefunctions, /++ = /* + f a l and 
/~~ = fa — fa- For large separation of excitons, the 
exponential factors prohibit both configurations from si- 
multaneously contributing to the wavefunction. Thus, a 
node can be approximated in the scattering wavefunction 
by simply requiring that U a (r) be zero for all r > R n . 
The error introduced by this approximation is of order 
exp(— 2aR n ), and is another limit on the use of small val- 
ues for R n . Since we only do calculations for low energy 
scattering, large R n , the lack of a well-defined exciton- 
exciton separation distance for short distances does not 
matter. 

This method for calculating scattering properties is 
very sensitive to the energy spectra {E a (Rn)}. To get 
accurate energies, we do not try to construct and opti- 
mize elaborate variational wavefunctions, but rather use 
DMC to project the energy from trial wavefunctions of 
the form given in Eq. (|l3|). The coefficients 7, c/, c g , 
and Ch are chosen to obey the cusp conditions on the 
wavefunction for small particle separations. The s-wave 
envelope functions U a (r) are taken as solutions to an em- 
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FIG. 2: Change in the projected energy eigenvalues 
A a (t) - A a (0) of the DMC states relative to the VMC 
states as a function of DMC projection time t, for basis 
states a = 1, 2, 3, 4. The eigenvalue equation is given by 
Eqs. (§)-(§), where the Hamiltonian and overlap matri- 
ces have been sampled using Eqs. ([ll]) and (p"2|). 

pirical exciton-exciton scattering potentials, 



V(r) = 



-Vo 1 - 



1/ d 



(14) 



■^4 • r <d 

r > d, 

where Vo and d have been self-consistently fit to ap- 
proximate the energy spectrum of the four particle scat- 
tering states. We take the guiding function ip to have 
the same form as the / ++ wavefunctions with U^(r) — 

[ce~ r2 / 8 + Eo^^qWI 1 '' 2 - The parameters {c,d a } are 
chosen to bias sampling towards the collision: {c = 
0, g?i = 2, di = c?3 = e?4 = 1} for ++ states and 
{c = 2, d a = 1} for (\>~~ states To check for convergence 
of the energies in DMC, we plot the energy difference, 
Ed mc ~ Ey mc ■ a s a function of projection time in Fig. |^. 
We see convergence after a projection time of 3 E" 1 . 

We thus find two energy spectra for each value of R n , 
as shown in Figs. 0(a) and ||(b). The spectra for the 
symmetric states ^ , show a clearly bound biexciton 
state, as seen in Fig. 0(a). The antisymmetric states 
4>~~, as shown in Fig. g(b), have no bound state. The 
binding energy of the biexciton is E-q = —0.0642(3) E x , 
in agreement with other ground state calculations, and 
is insensitive to the position of the node Rn because it 
is localized. In contrast, the delocalized scattering states 
are quite sensitive to i?jv, and their dependence on Rn 
is a measure of the elastic scattering matrix elements. 



D Calculation of phase shifts 

Using Eq. (||), we determine scattering phase shifts, 
5 ++ (k) and 5~~(k), which are shown in Fig. |] for the 
equal-mass case. The k — limits show us that there 
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FIG. 3: Energy spectra as a function of nodal position 
Rn, for (a) symmetric </> ++ states, and (b) antisymmetric 
<j> states, with m e — rrih- The lowest energy curve in 
(a) is the biexciton with binding energy E-q = 0.0642(3). 
These functions E(Rn) determine the phase shifts S ++ 
and S , by the relationship in Eq. (ph. 




FIG. 4: Phase shifts S(k) for the two s-wave scattering 
states, for m e = m h , calculate using Eq. (||) and the data 
from Fig. ||. 

is one bound symmetric state </> ++ and no bound anti- 
symmetric states </>~~. The slope at k = is related to 
the scattering length a s by 5'(k) = —a s . From a cubic 
polynomial fit to the data, with coefficients given in Ta- 
ble |, we find a = 4.57(2)0^ for the symmetric state, and 
a = 1.51(3) for the antisymmetric state. These model 
values are consistent with the measured upper bound Aa x 
found in CU2O. 



TABLE I: Coefficients for polynomomial fit to the low 
energy part of the phase shift functions for the case m e = 
m h . 
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S ++ tv -4.574(21) 


2.995(70) 


-0.829(54) 


8~~ -1.512(29) 


-0.138(90) 


0.216(64) 



TABLE II: Matrix elements for changing spin basis in the 
two exciton problem. Columns are in the |s, s ee , Shhjeh 
basis and rows are in the ex (s, s e h, s e /t| basis. 
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E Extracting spin dependence from spatial 
symmetries 

Single excitons can be in either a spin triplet or spin 
singlet energy eigenstates. (These states are degerate in 
our model, but are split by electron-hole exchange inter- 
action in more realistic models.) During collisions, how- 
ever, the spin of individual excitons is not a conserved 
quantity — the excitons can change their internal spin 
states during collision. To relate our scattering states 
4> ++ and <f> to cxcitonic spin states, we invoke the anti- 
symmetry of the total (spatial plus spin) fermionic wave- 
function. Fermionic antisymmetry implies that the states 
4> ++ and <\> have spin eigenstates |s00) e h and |sll) e ?i 
where the notation \ss e Sh) e h stands for total spin s, total 
electron spin s e and total hole spin s e . We denote the ex- 
perimentally relevant "excitonic spin basis" as \ssiS2)ex, 
where the total spin is s and the spins on the scattering 
excitons are s\ and S2, which take the values for singlet 
excitons and 1 for triplet excitons. Table [n] lists matrix 
elements for the change of basis \ss e Sh) e h —> \ssiS2}ex- 
The states |000) e / 1 and \sll) e h act as channels for s-wave 
exciton-exciton scattering. The matrix elements for scat- 
tering through these channels are: 



(k) = e /i(000; /c|s|000; k) e h 



1 



s—(k) = eh (Sll;k\s\Sll;k) eh = 



2ik 

^2i8 (k) 



1 



2ik 



(15) 



(16) 



The collision of two excitons with well defined initial 
and final spin states can be determined by decompos- 
ing the scattering events into the two channels, s ++ (fc) 
and s~~(k). Using the change of basis matrix (Table pi), 
we find s(k) = a + s ++ (k) + a_s (fc), where the co- 
efficients a + and a_ for all non -zero s-wave scattering 
processes are given in Table III. The s-wave scattering 
cross sections are given by a(k) = 87r|s(fc)| 2 , where there 
is a factor of two enhancement due to the identical parti- 
cle statistics. The spin-dependent cross sections take the 
form, 

a(k) — — [c++ sin 2 <5 ++ + c__ sin 2 



k 2 

-c+- 



sm 



'(*++-*-)], 



(17) 
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TABLE III: All non-zero s-wave scattering process. Ini- 
tial and final states are denoted by \ssiS2}ex, where s 
is the total spin and s\ and S2 are the individual ex- 
citon spins. The coefficients a+,a_ are for comput- 
ing the scattering matrix elements, and the coefficients 

c ++ , c ,c_| are for computing cross sections. The last 

column list the ratio of the scattering length a s to the 
cxciton radius a x for the case m e = m^- 
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where c ++ = a + (a + + a_), c__ = 

c_| = a + a_ depend on the initial and final spin states 

and are tabulated in Table H. The scattering lengths a s 
are given by a s = — a + 5 ++ (0) — a_<5~~ (0), where the 
derivatives of the phase shifts are determined from the 
linear coefficients in Table [j[ 



IV. RESULTS 

A Spin dependence of scattering cross sections for 
the equal mass case 

The calculated spin-dependent scattering lengths for 
the case m e = are presented in the last column of Ta- 



ble |I|. These are the low-energy limit the cross-sections 
shown in Fig. ||. In Fig. [| we have plotted the s-wave 
scattering cross sections versus scattering momentum for 
the case m e = mh and all non-zero spin configurations. 
Fig. ||(a) shows scattering of two singlet-excitons. Scat- 
tering of two triplet-excitons is shown in Fig. |](b), where 
the solid line represents the spin aligned 5 = 2 state, and 
the dashed line represents the 5 = state. The 5 = 
state scatters particularly strong because it has a large 
contribution from the s ++ channel, which is enhanced 
by the weakly bound bicxciton. Triplet-excitons in a rel- 
ative 5 = 1 state are spatially antisymmetric and thus 
have no s-wave scattering. We show s-wavc scattering of 
triplet-excitons from singlet-excitons in Fig. ||(c). This 
state has two distinguishable excitons, and can scatter 
by both s-wave and p-wave processes. As can be seen 
in Table III, the only contribution to the cross section 
is from the weaker s channel. The coefficient for s- 
wave scattering is particularly small because only half 
the scattering process is symmetric (s-wave), and there 
is an additional factor of one-half which cancels the iden- 
tical particle factor. 

There is also an triplet- to singlet-exciton conversion 
cross section given in Fig. ||(d). Although this is an in- 
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exciton-exciton scattering with m e = m/j, for the pro- 
cesses: (a) singlet-singlet — > singlet-singlet, (b) triplet- 
triplet — > triplet-triplet for total spin s = (dashed line) 
and s — 2 (solid line), (c) triplet-singlet — > triplet-singlet 
with s = l, and (d) triplet-triplet — > singlet-singlet and 
singlet-singlet — > triplet-triplet, both with s = 1. 

elastic process in experimental situations, it conserves 
energy according to our model Hamiltonian because we 
do not have an explicit interband spin coupling. The con- 
version of two triplet-excitons to two singlet-excitons can 
be understood as an inter-exciton exchange of a pair of 
electrons (or holes) . Since the spins of the individual exci- 
tons do not correspond to symmetries of the Hamiltonian, 
they need not remain constant during scattering. This 
conversion process is a physical consequence of the two 
inequivalent scattering channels s ++ aad s . This e 
feet has been reported in experimentaEa and theoretica 
work on exciton scattering in quantum wells. 
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B Mass dependence of scattering lengths 

The dependence of the cross sections on mass may be 
numerically studied by our methods. In Fig. ^ we show 
our calculated af + and aj~ as a function of mass ratio, 
mh/m e - We find that the scattering length is remark- 
ably insensitive to the mass ratio for m e /mh < 10 (cor- 
responding to a wide range of semiconductors), but then 
diverges near m e /mh = 15. This feature is lost in previ- 
ously published theoretical treatments of exciton-exciton 
scattering. The divergence in near mh/m e = 15 is 
due to the acquisition of a biexciton vibrational state. 
For m e = nih the biexciton has no bound excited states, 
while a H2 molecule has 15 bound vibrational states. 
Our calculations have shown the first of these appears 
near m^/m e = 15, with dramatic effects on the scat- 
tering length. The a~ s curve is relatively featureless 
because there are are no bound antisymmetric states in 
this range. We interpret the upward drift of aj~ for 
larger mass ratios as a systematic error due to difficul- 
ties in projecting states in excitonic systems with very 
different electron and hole masses. The heavy particle 



determines the projection time (<■ 



l R~ 2 ) while the 



light particle determines the diffusion time step r. The 
difficulty in handling large mass ratios make the method 
(as presented here) complementary to calculations that 
use the Born-Oppenheimer approximation. 

It is important to realize that similar relationships 
must exist between the scattering length and other mate- 
rial parameters, such as the Luttinger-Kohn parameters 
describing realistic hole states, external strain, and spin- 
orbit coupling, to name a few. Theoretical studies of such 
effects will need similar high-accuracy scattering calcula- 



tions, but applied to more accurate Hamiltonians, and 
are an area for future research. 



V. CONCLUSION 

To summarize, we have shown that there are several 
significant elastic scattering processes for excitons, and 
have given numerically exact values for a widely used the- 
oretical model. We find strong triplet-triplet and singlet- 
singlet scattering, with weaker triplet-singlet scattering 
and triplet-triplet to singlet-singlet conversion processes. 
Scattering is relatively insensitive to the mass ratio for 
m e /mh i5 10, but becomes very sensitive and actually 
diverges near m e /m^ ss 15. DMC has been found be a 
good tool for this four-particle excited state calculation, 
since the detection of weakly bound states requires very 
accurate evaluation of the correlation energy. 

This computational approach should be extended in 
many ways. The extension to higher angular momen- 
tum states would give important corrections at higher 
scattering energies. Application to biexciton-biexciton 
scattering are possible, but would be a bit more difficult 
because the scattering wavefunction would then have to 
describe eight interacting particles. Most importantly, 
the method should be adapted to better Hamiltonians so 
that the sensitivity of the scattering length to material 
properties for a wide range of materials can be studied. 
This approach should be quite useful for quantum well 
problems, which have similarities to the bulk problem 
study here, but have many more experimental parame- 
ters that can affect exciton-exciton interaction. 
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